function [ H ] = MockH5CA2B2C( t )

    [T, Y] = ode45(@rigidc, t, [1 0 0 0 0]);

    H = Y';

end

function dy = rigidc(t, y)
    dy = zeros(5,1);    % a column vector
    k1 = 1;
    k2 = 0.5;
    k3 = 1;
    k4 = 0.5;
    
    dy(1) = - k1*y(1);
    dy(2) = k1*y(1) - k2*y(2);
    dy(3) = k2*y(2) - k3*y(3);
    
    dy(4) = k3*y(3) - k4*y(4);
    dy(5) = k4*y(4);
end
